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Abstract. - In studying the end-to-end distribution function G(r, N) of a worm like chain by 
using the propagator method we have established that the combinatorial problem of counting 
the paths contributing to G(r, N) can be mapped onto the problem of random walks with 
constraints, which is closely related to the representation theory of the Temperley-Lieb algebra. 
By using this mapping we derive an exact expression of the Fourier-Laplace transform of the 
distribution function, G(k,p), as a matrix element of an inverse of an infinite rank matrix. 
Using this result we also derived a recursion relation permitting to compute G(k,p) directly. 
We present the results of the computation of G(k, N) and its moments. The moments < r 2n > 
of G(r,N) can be calculated exactly by calculating the (1,1) matrix element of 2n-th power of 
a truncated matrix of rank n + 1. 



The theory of flexible polymers is now understood [|lj- ||. Many polymer molecules have 
internal stiffness and cannot be modeled by the model of flexible molecules developed by 
Edwards. This is especially true for several important biopolymers such as actin, DNA, and 
microtubules ||. Models of semiflexible polymers have also applications in different topics 
besides polymer physics ||. If the chain length decreases, the chain stiffness becomes an 
important factor. A quantitative measure for the stiffness of the polymer is the persistence 
length l p , which is the correlation length for the tangent-tangent correlation function along the 
polymer. Polymers with contour length L much larger than l p are flexible and are described 
by using the tools of quantum mechanics and quantum field theory |l|- The standard 
coarse-graining model of a wormlike or a semiflexible polymer was proposed by Kratky and 
Porod [Q. A few first moments of G(r,N) were computed in ||- pJ. The literature on the 
computation of G(r, N) and its moments can be found in the book by Yamakawa [jllj. For 
recent work see |Q- Q. 

In this Letter we will study the problem of computation of the distribution function G(r, N) 
of the end-to-end distance of a semiflexible polymer, which is described by Kratky-Porod 
model, by using the analogy of the worm like chain with the quantum rigid rotator in an 
external homogeneous field || , |l5| within the quantum mechanical propagator method |l6| . 
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The relation to the rigid rotator was used in the recent work |14| . Relating the combina- 
torics of counting the paths contributing to G(r, N) to random walks with constraints and 
the Temperley-Lieb algebra [ p"7| , we have derived an exact formula for the Fourier-Laplace 
transform of the end-to-end distribution function, G(k,p), as an inverse of an infinite rank 
matrix of a very simple structure, given by Eqs.(^-0). Using a truncated matrix of rank n + 1 
one obtains n moments of r 2 exactly from the (1,1) element of the 2n-th power of this matrix, 
i.e., by applying this matrix 2n times to a canonical basis vector. Moreoever, this enables one 
to compute the distribution function, which describes n moments of r 2 exactly and all higher 
moments approximately. The truncation corresponds to taking into account the intermediate 
states with the quantum number of the angular momentum I < n (see below). 

The Fourier transform of the distribution function of the end-to-end polymer distance of 
the Kratky-Porod model [§ G(k, L) = / d 3 i?exp(-ik(R - R ))G(R- R , L) is given by the 
path integral as follows 

G(k,L) = J Dt(s)Y[d(t(s) 2 -l)eM-^J o L dst(s)-± j\(^-) 2 ), (1) 

where l p is the persistence length, and t(s) is the tangent vector at the point s along the contour 
of the polymer. The product over s in Eq.(Q) takes into account that the polymer chain is 
locally inextensible. If the polymer interacts with an external field, the potential energy 
— J Q L dsV(r(s)) should be added in the exponential of Eq.([j]). Notice that the representation 
of G(k, L) by using the cumulants of the moments of the end-to-end distribution function 
G(k, L) — exp(— J2m=i Jo c ^ s / i 2m(s)(k 2 ) m ) permits to write the following differential equation 
for G(R, L), d G(R, L)/dL ~ E^ =1 (-l) m+1 /i2m(^)AgG = 0, where A = V 2 is the Laplace 
operator. The latter with the term associated with the external potential was used in |L8| to 
describe the localization of a semiflcxible polymer onto interfaces and surfaces. 

The path integral (Q) (without the term depending on k) corresponds to the diffusion of a 
particle on unit sphere, |t(s)| = 1. This problem is equivalent to the Euclidean rigid quantum 
|. The Green's function of the rigid rotator or a particle on unit sphere fulfills 
the following equation 

^Po(0, V, L; (h po, 0) - ^l v Po = S(L)5(Q - O ), (2) 

where f2 is the spheric angle associated with angles 9, and if, and 5(£l— fio) is a two dimensional 
delta function. Henceforth instead of the contour length L we will use the number of statistical 
segments N = L/l p . Notice that in Quantum Mechanics N corresponds to the imaginary time 
it. The solution for Po(0, ip, N; 6o,<Po, 0) is 

P o (0, <p, N; O , Vo, 0) = E CX P(~ + 2 1)N ) Y ^(0, vmM, V0 ), (3) 
i,m 

where Yi m {0, <p) are the spherical harmonics. At a given I, m takes the values — I, —I + 1, 
I. The quantity Po(0, tp, N; 6q, tpo, 0) is given by Eq.([l]) with k = and with the following 
boundary conditions for the path t(s) (0 < s < N): t(N) = (0, <p), and t(0) = (9o,<po)- 

We now will consider the Green's function P(0, (p, N; 0o, <Po, 0) associated with Eq. ([!]). The 
differential equation for P is 

Jj^P(0, ¥>, N; 0o, <p , 0) - iv 2 ^P + U(Q)P = d(N)d(Q - Sl ), (4) 



rotator M, |15 
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where U(Q) = ikfi =ikcos(9), and fc is measured in units of l p . As it is well-known from 
the propagator method ]T^ | the latter can be written as an integral equation as follows 



N 



ds / dVl' P (O, N-,n' , s)U(n')P(n' , s;fi o ,0). (5) 



Eqs.([|||) describes also the Euclidean rigid quantum rotator in an external field. The it- 
eration of the latter generates the perturbation expansion of P in powers of the poten- 
tial [/(£!). The Fourier transform of the distribution function of the end-to-end distance 
of the polymer is obtained from P(fl, N; flo, 0) by integrating the latter over SI and ilo: 
G(k,N) = Jdnj dn o P(n,N;n o ,0). The factor 4tt is extracted from P(0, TV; fi , 0). Eq.(|) 
with fixed SI and Slo is the distribution function of the end-to end distance with fixed tan- 
gents of the polymer ends. The coefficient in front of (fc 2 )" of the expansion of G(k,p) 
in powers of fc 2 , multiplied by the factor (— l)™r(2n 4- 2) is the Laplace transform of the 
moment < r 2 " > of the end-to-end distribution function. Due to the folding character of 
the integrations over the contour lengths in (||), the Laplace transform of (||) with respect 
to N permits to get rid of integrations over the contour length. Similar to the propagator 
method used elsewhere |T(| each order of the perturbation expansion can be represented by 
graphs, examples of which are schematically shown in Fig.l. Each piece of the continuous 
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Fig. 1 - a) The 6th order of the perturbation expansion of the end-to-end distribution function, b) 
Random walk with constraint associated with the last line in a). 

line between two consecutive vertices are associated with the Laplace transform of the bare 
Green's function given by Eq.(||). The integral associated with each vertex is (cos#) ; , l , = 
J dQYy m ,(9, (p) cos 9Y^ m (9, tp). Using the property of the spherical functions cos 9Y^ m = 
y /((l + l) 2 - m 2 )/(4(l + l) 2 - m +1<m + y/(/ 2 - m 2 )/(4? 2 ~T)yi_ l! m gives {cos 6) v lm , m - 

y/((l + l) 2 - m 2 )/(4(Z + 1)2 - l)8 v , l+1 5 m ., m + V(' 2 - m 2 )/(4Z 2 - 1)^-!^,™. 

Due to the integrations over and f2 the quantum numbers I and m associated with the 
transition probabilities (amplitudes in the quantum mechanical counterpart of the problem) 
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at graph ends are equal to zero. The averages over cos0 in each vertex gives that the quantum 
numbers I on both sides of a vertex differ by one. The weights associated with the vertices are 
given by wi — (cos9) l l _ 1 = (cos8) l _ 1 1 — Z/a/ 4Z 2 — f . The Laplace transform of the weights 
associated with each piece of the continuous line between two vertices, which are given by 
the exponential factor in Eq.(^), are given by e; = + 1(1 + l)/2), where p is the Laplace 
conjugate to N. The quantum number m is conserved along the path and is zero due to 
the integration over the angles at the ends. Thus, the problem of counting different paths 
contributing to the transition probability J cZfi J drioP(fl, N; Qq, 0) can be formulated as a 
problem of a random walk consisting of odd number of steps with the constraint that the 
step heights are equal to one. The example of such a random walk is given in Fig. lb. The 
computation of the coefficients in front of (k 2 ) n of the expansion of G(k,p) in powers of k 2 is 
thus related to the enumeration of configurations of such a random walk consisting of 2n + 1 
steps. This kind of the problem is closely related to the calculation of the dimensions of the 
irreducible representations of the Temperley-Lieb algebra [jl7|] . One can write the expansion 
of G(k,p) as follows 



G(k,p) 



J2(-k 2 TA 2 



n=l 



[I + k 2 A 2 ] 



1,1 



(6) 



i,i 



where the last term is the result of the summation of the series, A = d.m is the product of 
the infinite rank matrices d and m, which are defined as follows d^j — ej_i#ij and rriij — 
Wi-i5i.j+i + Wi5i.j-i (i,j = 1, 2, ...). Approximating the last term in Eq.(||) by the matrices 
of finite rank n gives an approximate formula for the distribution function, C?2rt-2 (&>£>)) which 
gives exactly the first 2n — 2 moments and approximate expressions for higher moments where 
the intermediate states with I = 0, 1, ...,n — 1 are taken into account. No matrix inversion is 
required for the exact Laplace transforms of the moments which may be calculated directly by 
the Fourier transform of (||). From the expansion in k one realizes that only the (l,l)-element 
of the 2n-power of A needs to be calculated for < r 2n > . Because of the simple band structure 
of A this is obtained exactly from a truncated matrix A of rank n + 1 . 

The distribution function G2n(k,p) (2n is the number of the exact first moments contained 
in G2n{k,p)) can be written as 

n (V \ 1 l + /2 + /3 + - + /n ff7 , 

G 2n {k,p) = — — 7", (7) 

+ + J2 + h + ■■■ + f n 

where which are defined for each n, obey the following recursion relation fk = Wfe(l + 
-.k+2 fl)i an d u i — k 2 ei-\w 2 ei . Eqs.(|^-0) are exact and are the main result of the present 
Letter. The off diagonal matrix elements of the matrix [/ + fc 2 A 2 ] 1 give the transition 
probabilities from the initial state (l a , m a = 0) to the final state (lb, rrib — 0). The computation 
of the distribution function of the end-to-end distance with fixed tangents of the end vectors 
demands consideration of paths with arbitrary l a , m a and lb, vrib- The consideration of the 
transition probability with nonzero number m demands the change of the definition of the 
weights wi . The simple band structure of the matrix, however, remains unchanged. 

We have computed G m (k,p) explicitly for m < 26. The distribution function, which takes 
into account exact 6 moments, is given by 

5 39 k 2 + 11 k 2 p + 21 p 3 + 210 p 2 + 378 + 567 p 

^ ' P ' ~ 3 35 p 4 + 350 p 3 + 30 k 2 p 2 + 945 p 2 + 630p + 170 k 2 p + 3 k 4 + 210 k 2 ' 
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As an example we give below the 10th moment 
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All 26 moments we have computed are obtained in a form similar to < (r 2 ) 5 > with numerical 
coefficients given as fractional numbers. All moments fulfil the expected limiting behaviour 
for small and large TV. 

Performing the inverse Laplace transform of G(k,p) by using Maple enables one to plot 
G(k, N). In the limit of small TV the distribution function G m (r, TV) approaches the rod limit 
G rod {r, TV) = l/(47rTV 2 )<S(r - TV), which has the Fourier transform G rod (k, TV) = sin(kN)/kN. 
In order to illustrate the accuracy of the truncated dustribution function a comparison of 
G(k, TV) with the rod limit is shown in Fig. 2. The increase of m extends the interval of k of the 
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Fig. 2 - Computation of G(k, TV) for TV = 0.01: Continuous: G ro d(k,N), dots: Gio(fc, TV) 

agreement between G m (k,N) and G ro d{k, N). However, for large m the Maple computation 
of the inverse Laplace transform of G(k,p) becomes unstable for small k at very small TV. 
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The inverse Laplace transform of G(k,p)/p 2 multiplied with the factor 2/N gives the 
structure factor of the semiflexible polymer S{k, N). Performing the inverse Laplace transform 
by using Maple permits to get plots of the structure factor. 

In conclusion, we have studied the end-to-end distribution function of a worm like chain, 
the problem which is equivalent to the Euclidean quantum rigid rotator in a homogeneous 
external field, by using the propagator method. We have established that the combinatorial 
problem of counting the paths contributing to the end-to-end distribution function can be 
mapped to the problem of random walks with constraints, which is closely related to the 
Temperley-Lieb algebra. Using this mapping we derived an exact expression for the Fourier- 
Laplace transform of the end-to-end polymer distance as the matrix element of the inverse 
of an infinite rank matrix. Starting with this result we have also derived a recursion relation 
permitting to compute G{k,p) directly. The mapping to the problem of restricted random 
walks and the Temperley-Lieb algebra, which is not restricted to the computation of G(k, N), 
gives a new tool of treating the semiflexible polymers. Eqs.(^-0) give also the exact solution 
for the transition amplitude of a rigid quantum rotator in an external homogeneous imaginary 
field in Euclidean sector. 
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